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Abstract 

We  examine  the  shrinking-core  model  of  aluminum  combustion,  one  that  has  been  proposed 
for  sub-micron  drop  diameters.  In  this  model  a  core  of  liquid  aluminum  is  surrounded  by  alumina 
through  which  O  atoms  diffuse  to  the  aluminum  surface.  It  is  shown  that  the  volume  changes 
intrinsic  to  this  model  necessarily  lead  to  fracturing  of  the  alumina  and  the  creation  of  cracks  and 
voids.  A  simple  mathematical  model  is  described  which,  because  of  the  length  and  time  scales, 
is  quasi-steady,  permitting  analytical  solutions.  When  the  drop  is  sufficiently  small  this  leads 
to  a  simple  formula  for  the  burn  time  as  a  function  of  the  atmospheric  pressure,  temperature, 
and  oxygen  concentration,  which  is  tested  against  experimental  data.  By  introducing  a  fractal 
ingredient  into  the  description,  motivated  by  the  fracturing,  it  is  possible  to  generate  agreement 
with  experimental  data  on  the  variations  of  burn  time  with  drop  diameter,  a  d025-t  law  in 
contrast  with  the  familiar  d2-t  law  of  classical  fuel-drop  combustion.  For  larger  drop  diameters 
a  non-linear  differential  equation  for  the  burn  rate  is  derived  whose  integration  yields  the  burn 
time  as  a  function  of  drop  diameter  and  shows  a  transition  from  the  d0  25  law  to  a  d2  law. 

keywords',  nano-aluminum,  shrinking-core,  fracture,  fractal,  combustion,  burn-time/diameter 
laws 


1  Introduction 

The  combustion  of  aluminum  particles  is  an  important  subject,  if  only  because  such  particles  are 
a  common  ingredient  of  heterogeneous  rocket  propellants.  A  large  literature  exists  for  particles 
of  diameter  greater  than  100  /un  or  so  (see  the  review  of  Beckstead  [1]),  and  the  basic  structure 
of  the  combustion  field  is  that  of  the  classical  fuel  drop.  It  is  well  known  that  the  so-called  d2 
law  is  then  relevant,  a  linear  relation  between  the  square  of  the  diameter  and  the  burn  time. 

This  follows  directly  from  a  simple  dimensional  argument,  one  that  we  shall  repeat  here  as 
pedagogical  texts  frequently  muddy  the  issue. 

Almost  all  diffusion  flames  are  characterized  by  a  large  Damkohler  number,  so  that  the 
chemistry  is  fast,  and  does  not  define  a  relevant  time  scale.  Accordingly,  there  are  only  two 
relevant  lengths:  the  drop  diameter  d\  and  the  diffusion  length  pD /M  where  M  is  the  mass  flux 
per  unit  area,  proportional  to  (d/d t)d.  It  follows  that  d  oc  1/M,  whence  (d/d t)d2  =  const. 

More  recent  research  efforts  have  been  concerned  with  sub-micron  particles,  even  particles 
on  the  nano  scale,  and  important  experimental  work  has  been  carried  out.  One  immediate 
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observation  is  that  a  reduction  in  size  decreases  the  diffusion  time  scale  and  so  for  sufficiently 
small  particles  the  fast-chemistry  assumption  will  be  false,  and  the  d2  law  will  no  longer  prevail. 
This  has  been  confirmed  experimentally.  Modeling  efforts  have  been  reported  by  Zachariah 
[2,  3], who  considers  what  is  called  the  shrinking-core  model. 

The  shrinking-core  model  is  spherically  symmetric.  At  the  core  is  a  sphere  of  liquid  aluminum 
(aluminum  is  liquid  between  933  K  and  3100  K),  and  surrounding  this  sphere  is  a  shell  of 
solid  alumina  (AI2O3,  melting  point  2345  K).  When  immersed  in  a  hot  atmosphere  containing 
oxygen,  some  of  the  latter  dissociates  heterogeneously  at  the  oxide  surface  and  oxygen  atoms 
diffuse  through  the  shell  to  react  at  the  aluminum  surface.  This  surface  retreats  with  time  and 
the  traveling  reaction  front  is  a  transition  zone  between  the  aluminum  and  the  alumina.  The 
diffusive  flux  of  oxygen  to  the  front  is  comparable  to  the  mass  flux  through  it  as,  overall,  it 
sustains  the  reaction 


2 Al  +  3 O  ->  Al203 

one  in  which  53.96  mass  units  of  aluminum  are  consumed  for  48  mass  units  of  oxygen. 


(1) 


2  Some  Geometry  Issues 


Here  we  discuss  some  geometrical  issues  that  arise  because  of  the  material  conversion  in  the 
drop/particle.  These  arise  by  consideration  of  the  mass  of  aluminum  and  oxygen  that  goes  into 
the  formation  of  the  alumina.  As  supplied,  aluminum  particles  have  a  thin  coating  of  alumina 
(2-4  nm  thick)  and  we  neglect  this  in  our  discussion. 

Consider  a  sphere  (aluminum)  of  initial  diameter  di  which  is  converted  to  a  sphere  of  alumina 
of  diameter  d / .  Then  the  initial  drop  mass  is 


the  final  mass  is 


and  the  added  mass  of  oxygen  is 


Mass  conservation  then  requires  {pai 


53.96  6  iPal' 

2700  kg/m3,  paia  =  3950  kg/m3) 


=  (1.2916)3  =  1.0890. 

CLi 

Thus  one  might  expect  an  8.9%  increase  in  diameter. 

However,  let  us  examine  the  neighborhood  of  the  front,  which  we  suppose  is  located  at  a 
radius  s  at  time  t.  This  front  separates  aluminum  from  previously  generated  alumina,  and  is 
moving  with  speed  u  into  the  aluminum.  At  the  same  time  the  boundary  between  previously 
generated  alumina  and  newly  generated  alumina  is  moving  with  speed  v  in  the  opposite  direction; 
this  speed  arises  because  of  the  volume  demands  of  the  alumina.  Then  in  an  infinitesimal  time 
At,  and  considering  unit  area  of  front,  the  mass  of  aluminum  consumed  is 


the  mass  of  alumina  produced  is 


PaiuAt, 
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and  the  mass  of  consumed  oxygen  is 


Paia(u  +  v)At, 


Thus 


48 

53.96 


PalaUAt. 


-  =  0.2916. 
u 

In  other  words,  the  local  (linear)  argument  shows  that  the  overall  particle/drop  diameter  in¬ 
creases  by  a  factor  1.2916  over  the  complete  burn,  a  29.16%  increase  rather  than  an  8.9% 
increase. 

To  reconcile  these  facts  we  note  that  although  alumina  is  a  valuable  ceramic  with  many 
desirable  properties,  it  is  brittle  and  weak  in  shear  and  tension.  And  so,  although  some  defor¬ 
mation  can  be  sustained,  significant  fracturing  will  occur,  creating  voids.  If  we  assume  that  all 
of  the  adjustment  is  made  in  this  way,  we  can  calculate  the  effective  final  density  of  the  alumina 
sphere  by  assigning  its  final  radius,  viz. 


whence 


4  _  2700  I"  48 
(If  peff  .  +  53.96 


1.29163 


Peff  =  2368. 

The  history  of  the  density  change  can  be  described  by  considering  a  partial  burn.  Consider  a 
point  at  a  radius  r  (<  d/2)  fixed  in  the  laboratory  frame.  As  the  front  passes  through  this  point 
it  processes  a  shell  of  aluminum  to  create  a  shell  of  alumina  of  thickness  <5  located  between  r  and 
r  +  6.  At  some  later  time  when  the  front  is  at  a  radius  s  this  shell  has  been  displaced  outwards 
by  a  distance  0.2916(r  — s)  so  that  it  is  located  between  r+0.2916(r  —  s)  and  r  +  0.2916(r  —  s)  +  <5. 
Its  initial  volume  is  4irr25  and  its  final  volume  is 

47t  [r  +  0.2916(r  —  s)]26, 

so  that  the  final  density  is 

3950[1  +  0.2916(1  —  s/r)]~2. 

Figure  1  is  a  plot  of  the  density  ratio  pratio  =  p/3950  vs.  2 r/d  for  different  values  of  2 s/d. 
Consider,  for  example,  2 s/d  =  1/4.  At  2 r/d  =  1/4  the  density  ratio  is  equal  to  1,  and  as  r 
increases  it  falls  monotonically  to  the  value  0.6733  at  2 r/d  =  1,  the  drop/particle  boundary. 
The  physical  location  of  this  boundary  is  greater  than  d/2  of  course,  because  of  the  expansion. 
If  r'  is  the  radial  measure  once  expansion  is  started  a  point  at  r  <  s  is  mapped  onto  r'  =  r;  but 
a  point  at  r  >  s  is  mapped  onto  r'  =  r  +  0.29(j —  s). 

One  might  reasonably  ask  whether  a  ceramic  of  which  only  60%  is  solid  and  in  which  the 
void  space  is  created  by  cracking  and  displacement  can  have  any  structural  integrity,  would  not, 
in  fact,  disperse  in  fragments.  Yet  the  latter  does  not  occur,  [4].  One  possibility  is  that  the 
liquid  aluminum  is  drawn  into  the  cracks  by  capillary  action,  a  process  of  system  self-healing, 
although,  conceivably,  outer  layers  could  be  shed.  If  aluminum  is  drawn  up  into  a  crack,  gases 
must  be  pulled  down  through  another  crack  to  replace  the  volume.  And  when  the  burn  is 
complete  the  alumina  will  have  a  hollow  core.  Hollow  post-burn  cores  have  been  reported  by 
Zachariah  [3]. 
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Figure  1:  Density  ratio  praUo  vs  2 r/d  for  different  2 s/d. 


3  Mathematical  Model 

In  this  section  we  describe  the  mathematical  model  and,  from  scaling  considerations,  describe 
rational  approximations  to  the  equations.  A  key  observation  is  that  the  temperature  field  and 
the  oxygen  concentration  field  are  quasi-steady  so  that  the  problem  reduces  to  a  non-linear 
differential  equation  for  the  location  of  the  reacting  front  as  a  function  of  time.  If  we  take  some 
liberty  with  the  diffusion  coefficient  this  equation  reduces  to  a  simple  expression  for  the  front 
speed  in  terms  of  the  atmospheric  conditions  (temperature,  pressure,  oxygen  concentration). 
Certain  free  parameters  (e.g.  the  power-law  dependence  on  the  oxygen  concentration)  can  be 
chosen  to  fit  experimental  data,  and  for  a  fixed  particle  size  ( d  =  80  nm)  variations  of  burn  time 
with  the  oxygen  concentration,  and  burn  times  with  temperature  at  two  different  pressures  are 
accurately  calculated.  However,  the  predicted  burn  law  is  a  linear  one,  and  recently,  Glumac 
[5]  measured  the  burn  times  for  particles  of  size  18,  50,  80,  and  110  nm,  with  atmospheric 
conditions  fixed,  data  which  can  be  well  fitted  with  a  d 025  law.  And  so  we  have  created  an 
epistemological  model  by  a  simple  modification  of  the  derived  model  which  forces  the  d0'25  law. 
This  modification  is  partly  motivated  by  the  geometric  discussion  of  Section  2,  and  introduces 
no  new  parameters  (apart  from  the  0.25).  Predicted  burn  times  differ  by  negligible  amounts 
from  a  MatLab-generated  fit  to  the  experimental  data. 

Later  we  examine  the  solution  sans  liberties  with  the  diffusion  coefficient. 

3.1  Governing  Equations 

We  put  aside  the  geometrical  issues  of  Section  2  (specifically,  the  motion  of  the  oxide)  and  write 
down  simple  field  equations  for  the  temperature  and  the  concentration  of  oxygen.  There  is 
no  significant  convective  transport  in  the  laboratory  frame,  and  so  in  both  the  aluminum  and 
alumina  heat  conduction  is  governed  by 
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(2) 


hi  (vf ) 


Here,  T  is  the  temperature,  r  the  radial  coordinate,  t  time.  The  density  p,  the  thermal  conduc¬ 
tivity  A,  and  the  specific  heat  C  are  assigned  values  appropriate  to  the  medium.  Also,  in  the 
alumina  the  concentration  of  oxygen  (kg-mol/m3)  is  controlled  by 


dc  Id  /  2  <9c\ 
dt  r 2  dr  \  dr) 

where  D  is  the  diffusion  coefficient. 

At  the  front,  the  temperature  is  continuous 

Pit  =  0, 

where  (+)  refers  to  the  alumina  side  of  the  front. 

Energy  conservation  is  given  by 


dT 

dr 


=  QSM, 


(3) 

(4) 


(5) 


if  the  difference  in  the  specific  heats  is  neglected;  M  is  the  outward  mass  flux  through  the  front, 
and  Qs  is  the  heat  released.  And  the  diffusive  flux  of  oxygen  to  the  front  is  related  to  M  by 


W0D 


3  Wgl0 
Wal 


■M, 


a  consequence  of  the  stoichiometry  defined  by  equation  (1). 
If  the  front  is  located  at 


(6) 


r  =  s(t)  ( ds/dt  <  0)  (7) 

the  mass  flux  is  defined  by 

M  =  -pai  —  .  (8) 

The  front  moves  because  of  consumption  by  reaction  of  aluminum  at  the  surface  of  the  aluminum 
drop,  and  we  suppose  that  this  can  be  described  by  the  law 


ds 

dt 


T  \n 

— —  exp  (-E/RT  +  E/RTref), 

4  ref  J 


(9) 


where  p  is  the  pressure  and  the  subscript  ref  denotes  reference  quantities  to  be  specified. 
Finally  we  impose  boundary  conditions  at  the  outer  boundary  r  —  dj 2  of  the  alumina 


=  eCT(T4  _  T4  )  +  _  Too); 


dT 

dr 


(10) 


where  a  is  the  Stefan-Boltzmann  constant,  e  is  the  emissivity  of  the  alumina,  and  Nu  is  the 
convective  Nusselt  number. 
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Parameter 


Value 


Pal 

2700  kg/m3 

Pal  a 

3950  kg/m3 

c 

900  J/kg-K 

D 

10”7  m2/s 

Nu 

2 

Qs 

3.11  x  105  kJ/kg 

e 

0.41 

A al 

250  W /m-K 

A ala 

25  W/m-K 

O 

5.67  x  10-8  W/m2-K' 

Table  1:  Table  of  parameter  values. 


3.2  Scaling 

Here  we  examine  the  equations  and  come  to  conclusions  as  to  which  terms  are  important  and 
which  terms  are  not.  To  this  end  we  scale  using:  length  d/2;  time  r  =  burn  time;  speed  d/2r. 
Then  the  non-dimensional  temperature  equation  can  be  written 

dt  r2  dr  \  dr )  ’  ,  . 

4tA  ^  ’ 

“  ~  pCd 2 ' 

In  the  alumina,  with  r=  100/zs,  d= 80  nm,  and  other  data  from  Table  1,  we  find  that  a  is 
roughly  106  and  so  the  temperature  field  is  strongly  quasi-steady,  apart  from  any  (rapid)  initial 
transients.  Similar  conclusions  are  true  for  the  temperature  in  the  aluminum  (the  equivalent 
a  is  107)  and  for  the  oxygen  concentration  in  the  alumina  (the  equivalent  a  is  6  x  103).  Since 
the  r-dependent  solution  of  the  steady  equations  is  singular  at  the  origin,  it  follows  immediately 
that  the  temperature  in  the  aluminum  does  not  depend  on  r.  Moreover,  the  scalings  reveal  that 
the  temperature  gradient  on  the  oxidizer  side  of  the  front  is  essentially  zero,  as  is  the  gradient 
at  d/2,  and  so  the  temperature  in  the  oxidizer  is  likewise  independent  of  r.  The  situation  is 
not  so  simple  for  the  front  condition  (6)  involving  the  concentration  c,  for  there  we  find  a  rough 
balance  between  the  first  order  estimates  of  the  oxygen  flux  and  the  aluminum  flux. 


4  Solution 


As  noted  in  Section  3.2  we  may  not  assume  that  c  is  constant,  although  its  description  is  quasi¬ 
steady.  Thus 


c  c0  ^  B  / 1  2  \ 

Cref  Ore/  d  y  r  d  J 


(12) 


for  some  B ,  where  c0  is  the  value  of  c  at  the  boundary,  r  =  d/2.  At  s,  dc/dr  oc  B  oc  ds/dt  so 
that 

r,  31 VglaPal  2ds_  (T3l 

WalW0DcrefS  dt' 


Thus  c  at  r  =  s  is  determined  as  a  functional  of  s  and  the  regression  law  can  be  written 
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ds 

dt 


=  -K 


P 

Pref 


Q 


T 

Tref 


n 

exp 


E 

RT 


E 


Co 


RTref 


Cref 


(  s  _  2s^\  3W ala  Paid?  d  /  S\' 

\d  d2  )  AWaiW0Dcref  dt  \d) 

(14) 


We  do  not  have  a  great  deal  of  confidence  when  assigning  a  value  to  D.  Measurements  indi¬ 
cate  extraordinarily  large  variations  with  particle  size  and  temperature,  variations  that  are  not 
understood.  And  there  is  no  consensus  on  what  the  value  of  D  should  be,  for  example,  for  an 
80  nm  particle  at  1500  K.  The  difficulty  is  that  D  can  not  be  measured  directly,  but  only  by 
inference,  and  it  is  not  clear  what  other  ingredients  might  affect  the  measurements.  And  so  it 
would  be  unwise  to  confine  our  discussion  of  equation  (14)  to  the  specific  value  of  D  defined  in 
Table  1,  a  value  proposed  in  [6].  We  note  that  if  D  is  large  (10-5  -  10-6  say)  the  quantity  in 
brackets  can  be  approximated  by  c0/cref  and  we  have 


It  is  worth  noting  that  if  transport  of  oxygen  to  the  front  is  facilitated  by  the  presence  of  cracks 
(cf  Section  2),  a  large  effective  D  could  well  be  appropriate. 


E 


E 


RT  RT, 


ref  \ 


(15) 


4.1  The  large  D  approximation 

Here  we  examine  equation  (15),  and  two  things  need  attention:  the  value  to  be  assigned  to 
c0;  and  the  mismatch  between  the  linear  d  —  t  relation  ( d 1  law)  predicted  by  (15)  and  recent 
measurements  by  Glumac  [5]  that  support  a  d0  25  law  for  particles  between  18  and  110  nm  in 
diameter. 

Consider  c0  first.  The  molar  concentration  of  O2  in  the  atmosphere  is 

n/V  =  xp/RT ,  (16) 


where  n  is  the  number  of  moles  in  the  volume  V,  p  is  the  pressure,  xp  is  the  partial  pressure  of 
oxygen,  R  the  gas  constant,  and  T  the  temperature.  Thus  with  p  measured  in  atmospheres,  T 
in  degrees  Kelvin 

n/V  =  12.18  ( xp/T )  kg-mol/m3.  (17) 

Thus  for  a  pressure  of  8  atm,  a  mole  fraction  (a:)  of  1/2,  and  a  temperature  of  1500  K,  we  have 
n/V= 0.0325  kg-mol/m3. 

How  this  is  related  to  the  concentration  in  the  outer  regions  of  the  alumina  is  a  question  we 
are  not  equipped  to  answer  at  the  present  time.  A  fraction  of  the  oxygen  molecules  dissociate 
on  impact  with  the  alumina  surface,  and  so  there  is  a  flux  of  O  atoms  to  this  surface.  The 
calculation  of  this  flux  also  presents  a  challenge  that  we  are  not  able  to  meet,  particularly 
since  at  8  atm  and  1500  K  the  mean  free  path  ~20  nm  and  so  the  Knudsen  layer  must  be 
accommodated.  Some  of  this  flux  attaches  to  the  surface  and  thence  migrates  into  the  interior; 
we  do  not  know  what  fraction  does  so. 

To  bypass  these  difficulties,  we  assume  that  c0  is  proportional  to  the  concentration  defined 
by  (16).  Then  with  the  meanings  of  K,  q  and  n  reassigned,  (15)  can  be  replaced  by 


^  (£)'(£)' (si)’-- 

A  further  modification  can  accommodate  the  experimental  results  of  Glumac,  results  which 
can  be  forced  by  replacing  (18)  by 
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again  with  a  reassignment  of  the  meaning  of  K.  This  can  be  regarded  as  an  epistemological 
model  (of  the  pragmatic  school),  but  can  also  be  viewed  in  the  context  of  Section  2.  Let  us 
write  the  regression  law  (9)  in  the  form 


d_ 

dt 


—A'kszK 


Cq 

Cref 


Pref 


T 


T, 


ref 


exp 


E 

RT 


E 

RTref 


(20) 


Now  suppose  that  the  fracturing  that  occurs  on  multiple  scales  as  the  burn  progresses  creates 
a  boundary  surface  between  the  aluminum  and  the  alumina  that  is  fractal  in  nature.  This  will 
have  no  affect  on  the  volume,  the  left-side  of  (20),  but  on  the  right  the  surface  area  (~  s2  for  a 
smooth  surface)  will  then  be  proportional  to  s2+v  (0  <  v  <  1)  for  some  v.  The  choice  in  (19) 
corresponds  to  v  =  0.75. 


4.1.1  Comparison  with  experimental  data 

When  (19)  is  integrated  (s  from  d/2  to  0,  t  from  0  to  r)  we  get 


1  /  d  \025/  x  \~m  (  P  VV  T  \  ”ex  \E_ 

K'\dref)  \Xref )  \Pref )  \Tref  )  ^  RT 


E 

RTref 


K'  =  Kd 


-0.25 

ref 


2~ 


1.75 


(21) 

and  we  compare  this  formula  with  the  experimental  results  of  [7].  We  start  by  taking  dref= 80 
nm,  pref =32  atm,  xref= 0.5,  Tref= 1460  K  and  with  d  =  dref,  p  =  pref,  T  =  Tref  so  that 
t  =  (1/ K')(x / xref)~m  we  obtain  the  curve  shown  in  Figure  2  when  K'  =  0.011/xs-1  and 
m  =  0.402.  Figure  2  also  shows  data  from  [7]  with  ±15%  error  bars,  as  suggested  by  our 
laboratory  experience  in  this  field.  The  measured  temperature  for  which  this  data  was  obtained 
is  1500  K  but  for  reasons  that  we  describe  in  the  next  paragraph  we  adjusted  that  to  1460  K,  a 
modest  shift  that  does  no  violence  to  plausible  error  possibilities,  Now  with  p  =  prefi  d  =  dref, 
x  =  xref  so  that 


f  T  ) 

—n 

|  exp 

'  E 

E 

K'  1 

\  Tref  ) 

RT  ~ 

RTref  _ 

we  construct  the  lower  curve  of  Figure  3  with  the  choice  n  =  0.4,  E  =  56000  J/mol-K.  The  upper 
curve,  for  a  pressure  of  8  atm,  is  obtained  by  choosing  q  =  0.6.  Again,  the  experimental  data  is 
from  [7].  It  is  in  examining  these  two  sets  of  data  that  we  were  led  to  the  choice  Tref  =  1460 
K,  for  with  a  choice  of  1500  K  the  fits  are  inferior  unless  we  use  different  values  of  E  for  the 
different  pressures.  There  can  be  no  objection  to  such  a  strategy  for  there  is  no  reason  to  believe 
that  the  set  of  true  reactions  which  lead  to  an  effective  E  at  32  atm  lead  to  an  identical  effective 
E  at  8  atm.  But  we  are  inclined  to  use  a  single  E  strategy  if  only  to  avoid  the  problem  of  what 
value  of  E  to  choose  when  p=  20  atm,  for  which  the  data  of  Figure  4  are  relevant. 

Figure  4  shows  burn  times  vs  d  when  p  =20  atm,  T  =  1550  K,  and  x  =  0.2,  with  experimental 
data  from  [7].  It  is  this  data  that  led  us  to  force  the  d025  law. 


4.2  Solution  for  general  D 

In  general,  the  concentration  c  (equation  12)  can  not  be  approximated  by  ca  and  the  burn 
time  must  be  determined  by  integrating  equation  (14)  after  modification  to  allow  for  the  fractal 
boundary.  Thus  a  choice  must  be  made  for  cQ  and  so  we  start  with  further  discussion  of  that 
issue. 

It  was  noted  earlier  that  for  a  pressure  of  8  atm,  a  temperature  of  1500  K,  and  a  mole 
fraction  of  oxygen  of  0.5,  the  mole  concentration  of  O2  is  0.0325  kg-mol/m3.  It  follows  that 
the  number  concentration  is  0.196xl026  m-3.  Moreover,  the  mean  molecular  speed  ^756  m/s 
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corresponding  to  a  mean  speed  normal  to  the  surface  of  ~481  m/s.  The  flux  towards  the  surface 
is  therefore  ~7.8  kg-mol/m2s  and  so  the  flux  of  O  atoms  (albeit  bound)  is  15.6  kg-mol/m2s, 
equivalent  to  ~250  kg/m2s.  These  estimates  assume  that  the  molecular  motion  is  in  thermal 
equilibrium,  and  so  neglect  Knudsen-layer  effects,  but  provide  a  rough  estimate  of  the  maximum 
oxygen  atom  flux  available  for  reaction  with  the  aluminum. 

This  can  be  compared  with  an  estimate  from  equation  (6)  (the  right  side)  for  an  80  nm 
particle/drop  with  burn  time  100  /./s,  which  yields  3.06  kg/m2s,  two-orclers  of  magnitude  smaller. 
This  difference  between  the  gas-phase  flux  and  the  required  diffusive  flux  in  the  solid  reflects 
the  fact  that  only  a  small  fraction  of  the  oxygen  molecules  moving  towards  the  surface  end  up 
as  attached  O  atoms. 

Attachment  occurs  to  alumina  molecules  at  the  surface.  Since  the  mass  of  an  alumina 
molecule  is  16.89xl0-26  kg  and  the  alumina  density  is  3950  kg/m3,  the  molecular  number 
density  is  233.87xl026  m-3  so  that  the  molecule  spacing  ~0.35  nm.  Thus  the  number  of  atom 
receptor  sites  (i.e.  the  number  of  alumina  molecules)  on  the  surface  is  ~  8.16  x  1018  m-2.  These 
atoms  move  into  the  solid,  one  alumina  molecule  at  a  time,  and  to  estimate  the  concentration 
near  the  surface  we  assume  that  on  interior  sheets  parallel  to  the  surface  that  are  one  alumina 
molecule  apart  (i.e.  0.35  nm  apart)  there  are  ~  8.16  x  1016  atoms  per  square  meter.  Thus  the 
number  concentration  is  ~  2.33  x  1026  m-3,  corresponding  to  ~  0.39  kg-mol/m3.  We  use  0.4 
for  the  calculations  that  lead  to  figures  5  and  6. 

Consider  now  equation  (14)  with  p  =  pref ,  T  =  Tref  and  the  fractal  factor  added  to  K: 


—  =  - Ks 075 
dt 


Co 

Cref 


2s2  \  3 WglaPald2  d  /  s\ 

d2  )  AWaiW0Dcref  dt\d) 


(22) 


In  Section  4.1.1  we  solved 


ds 

dt 


=  —Ks 


0.75 


C’ref 


(23) 


with  cQ  — >  x  and  cre/  — >  xref,  xref  =  0.5  and  chose  K  and  m  to  match  experimental  data  for 
different  values  of  x.  x  was  a  surrogate  for  c0  -  we  assumed  that  c0  was  proportional  to  x.  We 
found  that  K'  =  0.011  /xs_1,  see  equation  (21).  The  same  choices  for  K  and  m  are  sensibly  used 
here,  therefore,  with  cre/  =  0.5.  Also,  we  take  co=0.4  kg-mol/m3. 

Figure  5  shows  the  solution  of  (22)  for  different  values  of  D.  Note  that  once  D  exceeds  10~7 
the  trajectories  are  essentially  indistinguishable.  Figure  6  shows  variations  in  the  burn-time 
with  particle/drop  diameter.  For  diameters  smaller  than  ~  300  nm  the  d025  law  is  satisfied; 
for  diameters  greater  than  ~  1  /jm  a  d2  law  is  satisfied.  The  d2  law  arises  by  setting  the  right 
side  of  (22)  to  zero,  corresponding  to  the  vanishing  of  c  at  the  front.  This  arises  when  d  is  large 
and/or  D  is  small  so  that  the  fast-chemistry  limit  is  achieved.  It  should  not  be  confused  with 
the  d2  law  of  classical  fuel-drop  burning,  as  the  physics  that  we  embrace  here  is  quite  different. 


5  Conclusions 

In  this  paper  we  have  examined  the  shrinking-core  model  of  aluminum  combustion,  proposed 
for  small  drops,  and  provided  strong  evidence  that  there  is  a  parameter  domain  for  which  the 
burning  rate  is  controlled  by  the  reaction  rate  at  the  aluminum  surface.  Specifically,  measured 
burning  rate  variations  vs  pressure  variations,  for  example,  can  be  correlated  with  pressure 
variations  of  the  reaction  rate.  We  have  suggested  that  the  alumina  generated  during  the  burning 
process  will  fracture  to  create  cracks  and  voids,  and  that  this  will  generate  a  fractal  boundary 
which  is  responsible  for  the  d0  25  burning  law  observed  for  sufficiently  small  drops.  For  larger 
drops,  but  ones  still  small  enough  to  fall  within  the  shrinking-core  framework,  we  have  described 
a  transition  between  the  d025  law  and  a  d2  law  as  the  drop  size  is  increased.  It  is  important 
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to  note,  however,  that  this  should  not  be  confused  with  the  classical  d 2  law  observed  for  large 
aluminum  drops,  as  the  physics  is  quite  different,  and  the  temperatures  are  much  lower.  In  the 
discussion  of  the  model  we  have  examined  the  flux  of  oxygen  atoms  in  the  gas  towards  the  drop, 
and  noted  that  this  flux  is  two  orders  of  magnitude  greater  than  that  needed  in  the  combustion. 
This  suggests  that  the  oxygen  atoms  at  the  surface  of  the  alumina,  which  subsequently  diffuse 
into  the  interior  towards  the  aluminum  surface,  only  occupy  one  hundredth  or  so  of  the  available 
surface  receptors,  the  alumina  molecules.  With  this  assumption,  and  assuming  that  the  in-depth 
spacing  of  the  diffusing  atoms  near  the  surface  is  equal  to  the  diameter  of  an  alumina  molecule,  a 
concentration  of  atoms  can  be  estimated  that  yields  plausible  results  when  used  in  the  modeling 
of  the  larger  drops. 
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Figure  2:  Burn  times  vs  oxygen  concentration,  model  (curve)  and  data  (symbols)  from  Figure  6  of 
[7].  p  Prefi  T  Ae/>  d  dref. 
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Figure  3:  Burn  times  vs  temperature  at  two  different  pressures,  model  (curve)  and  data  at  p  =  32 
atm  (square)  and  p  =  8  atm  (triangle)  from  [7],  Figure  5.  x  =  xref, 


Figure  4:  Burn  times  vs  particle  diameter,  model  (curve)  and  data  (symbols)  from  [5] ;  p  =  20  atm, 
T  =  1550K,  x  =  0.2;  the  dashed  curve  is  a  Matlab  fit  to  the  experimental  data,  d0'2586. 
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Figure  5:  Variations  of  s  with  time  for  various  D ;  from  right  to  left  D  =  10  10,  9 
P  Prefi  T  Tref- 


d  (nm) 


Figure  6:  Variations  of  burn  time  with  particle  diameter  d\  D  =  10  'm2/s,  cQ  =  0.4, 
T  =  Tref.  The  transition  is  from  a  d0'25  law  to  a  d2  law. 
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